Universal zero-frequency Raman slope in a d-wave superconductor 
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It is known that for an unconventional superconductor with nodes in the gap, the in-plane 
microwave or dc conductivity saturates at low temperatures to a universal value independent of 
the impurity concentration. We demonstrate that a similar feature can be accessed using channel- 
dependent Raman scattering. It is found that, for a d^2_y2-wa.ve superconductor, the slope of 
low-temperature Raman intensity at zero frequency is universal in the Aig and -B29 channels, but 
not in the Big channel. Moreover, as opposed to the microwave conductivity, universal Raman 
slopes are sensitive not only to the existence of a node, but also to different pairing states and 
should allow one to distinguish between such pairing states. 



PACS numbers: 78.30.-j, 74.62.Dh, 74.25.Gz 

The effect of impurity scattering has been of increasing interests in studies of High- Tc superconductors . For 
superconductors with an order parameter which has nodes on the Fermi surface, it is well known that impurity 
scattering can lead to a finite density of quasiparticle states at zero energy (gapless excitations). In particular, in 
the strong resonant scattering limit, quasiparticle states can be strongly localized as a result of the short coherence 
length and mean free path provided that the impurity concentration is low. This has some interesting observable 
consequences. Of equal interest is that, as first predicted by Lee while the effective impurity scattering rate 7 is 
quite different for different impurity concentrations and different scattering limits (for example, 7 Aoe"^""" in the 
Born limit and 7 ~ Ao(Aor)~2 in the unitary limit with l/r the normal-state scattering rate and Aq the maximum 
of the gap), the microwave conductivity saturates at low temperatures (cto = J^e^/TrmAo) and is independent of 7 (or 
the impurity concentration) . The experimental verifications of this universal feature gives unambiguous evidence that 
the order parameter in High-Tc materials exhibits nodes on the Fermi surface. 

In the context of thermal conductivity, Graf et al. have found a similar universality related to the microwave 
conductivity via the Wiedemann-Franz law. This has been confirmed in an experiment by Taillefer et al. |^] which 
measures the in-plane low-temperature thermal conductivity of YBa2Cu306.9 at different Zn substitutions for Cu. In 
the present letter, we demonstrate how one can also study these universal features by doing channel-dependent Raman 
scattering experiments. It is found, in a d^2_y2--wa,Ye superconductor, that the low-temperature slope of the Raman 
intensity at zero frequency arc universal in the Aig and B2g channels, but strongly dependent on the scattering rate 
(^ 7^) in the Big channel. As opposed to the microwave or thermal conductivity for which an entire Fermi surface 
average is taken and thus the universal feature is general for all the pairing states so long as gap nodes cross the Fermi 
surface (except for a scale factor), which channels saturate and are universal and which do not saturate in Raman 
scattering change with the gap symmetry. This channel-dependent universality exhibited in the Raman scattering 
thus gives one more method to test the symmetry of the pairing state in the High- Tc superconductors and is more 
powerful than microwave or thermal conductivity because of the additional selectivity involved in the Raman geometry 
which allows some information on the position of the nodes in the Brillouin zone to be obtained. This phenomenon 
should be of particular interest in the heavy fermion superconductors where the pairing states are considered to be 
more diverse. 

The saturation of the low-temperature microwave conductivity in a c?-wave superconductor, is a result of a can- 
cellation between the value of the impurity-induced density of states at zero energy and the quasiparticle relaxation 
lifetime and arises only if there exists nodes in the order parameter on the Fermi surface. The saturation in the slope 
of the low-temperature zero-frequency Raman intensity can be understood in a similar manner, with an additional 
channel dependent feature unique to Raman which arises from the different dependence on the chosen Raman ge- 
ometry vertices which pick up different contributions around the Fermi surface. In a dx2_y2--wave superconductor, 
the i32g-channel Raman vertex selects preferentially states along the diagonals of the Brillouin zone (where the nodal 
lines are) and hence probes directly the low-lying quasiparticle excitations. Consequently the net result is similar to 
what is seen in the microwave conductivity case. In contrast in the Big channel, the Raman vertex has maximum 
weight along the or ky axes and zero weight along the diagonals, consequently one is effectively measuring a finite 
gap and no universality is observed. 

The Raman intensity is proportional to the imaginary part of the zero-momentum limit generalized density response 
function 
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Xfriii'n) = -T'^ Tr[7(k)G(k, zw„ + w„))f(k, w„)G'(k, iw„)], (1) 

where Tr denotes a trace and 7(k) = T37(k) is the bare Raman vertex with fs the Pauh matrix and 7(k) = p'Ss • 

1*^ • Gj (effective mass approximation |^ ) . Here is the electronic dispersion relation of the superconducting layer 

and ei (e^) correspond to the polarizations of incident (scattered) photons. The renormalized Raman vertex T in (Q) 
is given by 

f(k,it'„) = f37(k) - faWcT ^ Tr[f3G(k', ia;„ + i;/„)f(k', w„)G(k', ia;„)], (2) 

where Vc is the Coulomb interaction. In Eq. (H), we have ignored the contribution to the vertex corrections due 
to the impurity potentials and two-particle pairing interactions and have included only the effect of the Coulomb 
interaction. For isotropic impurity scattering, it is sufficient to use the bubble diagram at small q, while the inclusion 
of the pairing interaction vertex correction is shown to have little effect on the Raman spectra and is particularly 
negligible at the low frequencies of interest. However, the effect of impurities is fully included in the single-particle 
Green's function G in Eqs. @) and (|). Substituting Eq. (|) into (0), one obtains 

XirilVn) = X-y-iVI-Vn) — (3) 

Xii(«i'«j - Vc 

where X7i(*^ri) is defined in the same way as X'^TiiVn) with F replaced by fa in (^, and so on. 

In terms of the particle-hole space, the single-particle Green's function is given by G~^(k, itUn) = ^'^nTo — ^MTz~^'kfl^ 
where d;„, ^k, and Ak are the impurity-renormalized Matsubara frequencies, electron energy spectrum, and gap. 
G is related to the noninteracting Green's function G(^^(k, iwn) = ia;„fo — ^kfa — Akfi via the Dyson's equation 
G~ ^ (k, iujn ) = "'^ (k, iun ) — 5^ (k, iw„ ) . We shall solve the self-energy S due to the impurity scattering. By expanding 
Ti{iujn) = "^aii^njTa (« — 0, 1, 3), One finds ja)„ — iojn ~ So, |k = Ck + S3, and Ak = Ak + Si. Employing the 

usual T-matrix approximation, the self-energy is then given by S(k, ia;„) = nir(k, k, ia;„), where rii is the impurity 
density and 

f (k, k', = z;,(k, k')f3 + ^;,(k, k")f3G(k", zo;„)f (k", k', za;„). (4) 

k" 

Here Wi(k, k') = (k'|ui|k} is the impurity potential. If we consider only isotropic impurity scattering [wi(k, k') = Wi], 
the T-matrix in (^) is left only with frequency dependence and can be solved to get T(icij„) = [1 — ViT3G{iujn,)]~^ViT'i 
with the integrated Green's function G{iuJn) = J2k ^(k, i^^n)- One can expand G{iujn) ~ Ga{ioJn)Ta (a = 0, 1, 3) 

with Ga{ii^n) = 1/2 X^k '^^[''"^e*' For a superconductor with particle-hole symmetry and an odd-parity gap 
which is the case for the d^2_y2 state, Gi(ia;„) = G^^icUn) = 0. This immediately gives 
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where c = 1/vi. The result of Si = is a reflection of the well-known result that in an unconventional superconductor 
with nodes in the gap and zero average, the gap remains unrenormalized due to the impurity scattering (Ak = Ak). 
Furthermore the effect of S3 is absorbed into the chemical potential as usual and consequently ^k = ^k- Equation (|^) 
is convenient for discussing both the weak (Born) scattering (c ^ 1) and strong (resonant) scattering (c <C 1) limit. In 
the normal state (Ak — 0), one can easily work out that iujn ~ ii^n + «(l/2r), where in the Born limit, the (isotropic) 
scattering rate l/2r = 2'KN{0)nivf , while in the resonant limit, l/2r = ni/2'KN{0). Here A^(0) m/^-Kh? is the 
density of states per spin on the Fermi surface. 

One can use a general Raman vertex 7s (k) = 7° -I- 7^/^(0) to classify different symmetry channels [| 10| denoted 



by s. Here 7° represents the isotropic and 7] represents the anisotropic part of 7^ and </> is the azimuthal angle in 
the x-y plane. In the case of a cylindrical Fermi surface, fs{4>) = cos{2(f)) for the Big channel, fs{<j>) = sm{2(f>) for 
the B2g channel, and fs{<l>) = cos(40) for the Aig channel. With the above Raman vertices and in the case of perfect 
screening {v^^ ^0), Eq (m is reduced to 
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where we have defined 
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In (^), the isotropic term (7^) is dropped as it cancels due to Coulomb screening. In the following we shall limit 
ourselves to a (ia;2_j£-wave superconductor with gap Ak = Aq cos(2(/)). In both the cases of the Big and B2g channels, 
the second term in Q vanishes since xj(w„) = 0. Technically this is due to the Fermi surface average Ak|^) — 0. 

However this condition doesn't hold in the Aig channel where the squared gap function has a component identical to 
fAi {(f>) — 4cos(4(/)) and hence (cos(40)|Akp) 7^ 0. As a consequence, the Coulomb screening only has an effect on 
the Aig channel intensity and has no effect on the Big and B2g channels. 

Based on the cylindrical Fermi surface approach, the second term of Eq. (^) in the Aig channel is shown to be of the 
same order as the first term in the zero-frequency limit In general, however, the effect of the second term is quite 
sensitive to the underlying quasiparticle energy dispersion relation and the issue regarding the effect of screening in 
the Aig channel remains an issue of considerable debate . On the other hand, experimental data seems to suggest 
that the first term of Eq. (^) can account well for the Raman intensity in the low-frequency regime. We thus drop 
the second term of (^) in our calculations of the slopes of the channel dependent Raman intensity at zero frequency 
in all channels. These are defined by (im^ D, + iS) 



dn 



4rm 



n 



(8) 



where the double prime denotes an imaginary part and the equal sign arises because X^'rl^ = 0) = 0. Using the 
spectral representation for the imaginary frequency Green's function, analytically continuing to real frequency from 
imaginary frequency {itUn to -\- iS) , and then performing the frequency sum and momentum sum (replaced by an 
integration Ek = 2iV(0) /^^ d( ^) gives 
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where uj± = iuni^ ^ ^^); = *<^n(w — ft ±id) and ^± = sgn^uj) ^ Coj. — A^; = sgn(ci; — fl) 
are chosen to have branch cuts such that Im^_|_, Im^^ > and Im^_, Imi^^ < 0. The index 1 in the vertex is dropped 
(7I ^ 7s) for simplicity. 

It is useful to compare Eq. (|^) for the zero-frequency Raman slope with a similar expression for the microwave 
conductivity (see, for example, Eq. (2) of Ref. |^). One finds that they are the same apart from an overall constant 
factor which appears in front of the expression and from a different angular function [in Raman scattering, the angular 
function fs{(j)) = cos(40), cos(2(/)), or sin(40) for Aig, Big, and B2g channels, while in the conductivity the angular 
function is usually Px = cos(0) for calculating CTq^ or py = sin(0) for calculating CTq^]. Also the appearance of term 
A^ in the second line of (||) is unique to Raman and occurs due to the different type of vertex (which is coupled to 
T3 in Raman and tq in the conductivity). The term proportional to A^, however, will drop out due to a cancellation 
between the two terms at zero temperature and contributes only a small amount at finite temperatures. 

We consider first the T = limit which gives [f{oj) — f {ui — ft)] / il « df{u!)/duj w —S{lu) when SI 0. This means 
that the uj integration in (^ is sharply peaked around the small uj region centered at a; = 0. Consequently we have 



S = A^(0)7.' 



(7^ + A^)t 
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where (• • •) denotes an average over the Fermi surface. This expression agrees with one obtained before by Devereaux 
and Kampf | p^ . Here the impurity- induced scattering rate in the superconducting state at zero- frequency is 7 = 
—iCj{ijj = 0) = iEo(^ = 0). The self-consistent results for 7 in the two different scattering limits Born and resonant 
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were solved for by Lee as mentioned earlier. Assuming that 7 <^ Aq (which requires that the impurity concentration 
rii be small in the resonant limit), we find at T = 



s = B2g or Aig 

(11) 
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As shown clearly in (|Tl|), the zero- frequency Raman slopes in both B2g and Aig-channel exhibit a universal saturated 
value at T = which is independent of 7, i.e., of impurity concentration - a feature first discovered by Lee in 
the context of the microwave conductivity. The reason the B2g and Aig channels share the same limiting value is 
simply because the square of angular functions [/ssg ('Z')]^ — sin^(20) = 1 — cos^(20) and [fAig{4')]'^ = cos^(40) = 
[1 — 2 cos^ (Sf/))]^ and the contribution due to whatever terms couple to cos(20) is small. In contrast in the Big channel, 
the zero-frequency Raman slope is proportional to 7^ (up to a logarithmic correction) and hence is strongly dependent 

on the impurity concentration. In the unitary limit, 7 t~^/^ ~ ^l^^i therefore the Big slope S ~ n^. We recall 
that for a system with tight-binding bands, the Raman vertex strength 7Aig and 751^ in ( |ll[) is proportional to the 
nearest- neighbor hopping, while 753^ is proportional to the next nearest-neighbor hopping ]14|]. 

If we take the ratio between the Raman intensities from superconducting and normal states, we find in the limits 
of 57 = and T = 0, x's/Xn ~ 7^/^o''' ^ig channel and ~ I/Aqt in the B2g channel. These differ from the 

expression given by Devereaux and Kampf [|l3| in that they use 7 for 1/r in the normal state. 

The finite but low-temperature (T < 7) limit is obtained by expanding ui± = ±1(7 -I- bui^) + aui at small w in (|^), 
where 7, a and h are constants and are found to be a ~ 1/2 and h ~ 1/(87) in the resonant scattering limit of primary 
interest here. Expanding the integrand in (^ to second order in lv leads to the finite-T result 



m72 / TT^T^, 
S - -dV 1 + ^7^ s^B2gO.Aig 
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While the universality is channel dependent, the variation ^ is found in all three channels. Equivalent results were 
given by Hirschfeld et al. ^ for the microwave conductivity and by Graf et al. for the thermal conductivity. 

As mentioned before, the channel-dependent universal zero-frequency Raman slopes are sensitive to differences in 
pairing states. For completeness, we consider some other pairing states of interest and focus only on the Big and 
B2g channels at T = 0. The opposite result to (|ll|) with Big and -823 channels switched, is obtained if the gap has 
d^y symmetry with = Aosin(2</)). In this case, one finds a universal feature in the _Big-channel Raman slope, 
but not in the B2g channel. This is because one sees mainly regions of maximum gaps in the i32g-channel, while in 
the i?ig-channel, one sees mainly the regions of zero gap. For a system with extended s-wave pairing state, i.e., with 
Ak = Aocos(4</)), we find that both Big and i?2g channels have the same universal Raman slope S — m^l Aq. 
In this case, the gap nodes appear at the angles of 7r/8, Stt/S and the equivalent and consequently both Big and i?2g 
channels see the same effective contribution from the gap node regions. 

Finally, we consider the case of a mixing gap with a d2.2_y2-wave component and a small isotropic s-wave component 
of weight a with the gap given by Ak = Ao[cos(20)-|-a], where a < 1. Physically this is somewhat similar to a system 
with orthorhombic band structure and a pure c?2;2_j,2-wave gap. In contrast to the pure dj.2_y2-wsLve case, the gap 
nodes arc now shifted away from the diagonals and consequently the Big Raman channel will have some contribution 
from the node regions. This means that the universal zero-frequency Raman slope feature will also be present in the 
Big channel and its weight will depend on how far the gap node have shifted off the diagonals, i.e., how much node 
contribution this channel picks up. To leading order, we find 5* = m'j'^/Tr^h^Ao in the B2g channel (same as the pure 
dx2-y2-wave case), while in the Big channel, 

= + ^'/^o ln(Ao/7)] . (13) 

TT^h Aq 

Therefore when a ^ 7/A0, the Big Raman slope has universality which however breaks down when a ^ 7/A0. In the 
unitary limit, 7 ^ and thus Ssig ^ m^"^ /n'^h'^ AQ{a'^ + crii). A study of the impurity concentration dependence 
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of Ssig may be used to extract the size of which in turn enables one to test how orthorhombic the system is. We 
note that in the d^^ -y2 +s-'wave case, (Ak) ^ and Si in Eq. (|) is nonzero and must be retained in the calculation. 
These complications are accounted for (see Ref. in obtaining the above result. Moreover, in this case, the second 
(screening) term of Eq. (^) contributes in the Big channel (but not in the B2g channel), though this contribution is 
expected to be small when the s-wave component is small. 

In conclusion, we have found that channel-dependent Raman scattering can be used to test a universal low- 
temperature behaviour in disorder high- Tc superconductors which in turn can reveal information on the pairing 
state in these materials. We studied the slopes of low-temperature Raman intensity at zero frequency in various chan- 
nel. Similar to what was found in the microwave or dc conductivity 0| and thermal conductivity for a (ij,2_j^2-wave 
superconductor with nodal order parameter along the diagonals on the Fermi surface, the universal feature holds in 
both Aig- and i32g-channel Raman slopes. Such a feature is not found, however, in the i3ig-channel Raman slope. 
Moreover, in contrast to microwave or thermal conductivity where the universal features are general to all the pairing 
states with nodes and certain reflection symmetry, the channel-dependent universality or lack thereof in Raman scat- 
tering is sensitive to difference in pairing states and hence allows one to further clarify the symmetry of the pairing 
states. In addition for a mixed order parameter with d^^-y^ and s symmetries, we have shown that Raman is able 
to give the amount of s-wave component which is not the case using microwave or thermal conductivity. Finally, 
we remark that, as recently shown by Branch [p^ , the strong inelastic scattering effect in the high- Tc materials may 
dominate over the elastic impurity scattering but will not destroy the universality in Raman slopes predicted here 
because this is a zero-temperature effect and at T = the inelastic scattering time becomes infinite. 
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